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The perturbative approach to the description of long wavelength excitations at high temperature 
breaks down near the critical point of a second order phase transition. We study the dynamics of 
these excitations in a relativistic scalar field theory at and near the critical point via a renormaliza- 
tion group approach at high temperature and an e expansion in d = 5 — e space-time dimensions. 
The long wavelength physics is determined by a non-trivial fixed point of the renormalization group. 
, At the critical point we find that the dispersion relation and width of quasiparticles of momentum p 

■ is tjp ~ p^ and Fp ~ — 1) uip respectively, the group velocity of quasiparticles Vg ~ p^^^ vanishes 

' in the long wavelength limit at the critical point. Away from the critical point for T > Tc we find 

(N , 

(^p ^ S, ^ [l + (pf)^^] ^ ^iid Tp ~ (z — 1) i^p with f the finite temperature correlation length 

^ , (, (X \T — Tc\~^ ■ The new dynamical exponent z results from anisotropic renormalization in the spa- 

^ ^ ■ tial and time directions. For a theory with 0{N) symmetry we find z — 1 + e f^^^^yi + C'(e'^).This 

dynamical critical exponent describes a new universality class for dynamical critical phenomena in 
quantum field theory. Critical slowing down, i.e, a vanishing width in the long-wavelength limit, 
' and the validity of the quasiparticle picture emerge naturally from this analysis. 
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INTRODUCTION 



The experimental possibility of studying the phase transitions of QCD via ultrarelativistic heavy ion collisions with 
the current effort at RHIC and the forthcoming program at LHC motivates a theoretical effort to understand the 
^ • dynamical aspects of phase transitions at high temperature. QCD is conjectured to feature two phase transitions, 
I the confinement-deconfinement (or hadronization) and the chiral phase transitions. Detailed lattice studiesj^] seem 
to predict that both transitions occur at about the same temperature Tc ^ 170Mev. 

While lattice gauge theories furnish a non-perturbative tool to study the thermodynamic equilibrium aspects of the 
transition the dynamical aspects cannot be accessed with this approach. 

In a condensed matter experiment the temperature is typically a control parameter and it can be varied sufficiently 
slow so as to ensure that a phase transition occurs in local thermodynamic equilibrium. In an ultrarelativistic heavy 
ion collision the current theoretical understanding suggests that a thermalized quark-gluon plasma may be formed 
I at a time scale of order Ifm/c with a temperature larger than the critical. This quark-gluon plasma then expands 
hydrodynamically and cools almost adiabatically, the temperature falling off as a power of time T{t) ~ Ti{ti/t)^^^ 
until the transition temperature is reached at a time scale '--^ 10 — 50fm/c depending on the initial temperature |2|. 

Whether the phase transition occurs in local thermodynamic equilibrium or not depends on the ratio of the cooling 
time scale tcooi ~ T(t)/T{t) to the relaxation or thermalization time scale of a fluctuation of a given wavelength p^^, 
ireiip)- If tcooi 3> trei then the fluctuation relaxes in time scales much shorter than that of the temperature variation 
and reaches local thermodynamic equilibrium. If, on the other hand, tcooi <C trei the fluctuation does not have time 
to relax to local thermodynamic equilibrium and freezes out. For these fluctuations the phase transition occurs very 
fast and out of equilibrium. Thus an important dynamical aspect is to understand the relaxation time scales for 
fluctuations. 

A large body of theoretical, experimental and numerical work in condensed matter physics reveal that while typ- 
ically short wavelength {p ^ T) fluctuations reach local thermal equilibrium, near a critical point long wavelength 
fluctuations relax very slowly, and undergo critical slowing down^, ^. A phenomenological description of the dynam- 
ics near a phase transition typically hinges on the time-dependent Landau- Ginzburg equation which is generalized 
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to include conservation laws|^, ^. In the simplest case of a non-conserved order parameter, such as in a scalar field 
theory with discrete (Ising-like) symmetry, the time dependent Landau Ginzburg equation is purely dissipative. 

While phenomenological, this approach has proved very succesful in a variety of experimental situations and is likely 
to provide a suitable description of the dynamics for macroscopic, coarse-grained systems such as binary mixtures, 
etc.|^. The phenomenological approach based on the time dependent Landau- Ginzburg equations which are first order 
in time derivatives seem to provide a suitable description of coarse-grained macroscopic dynamics in non-relativistic 
systems. However it is clear that such approach is not justified in a relativistic quantum field theory, since the 
underlying equations of motion are second order and time reversal invariant. 

In particular in the case of purely dissipative time dependent Landau Ginsburg (phenomenological) description^, Q, 
frequency and momenta enter with different powers in the propagators, and at the mean field (or tree) level, this 
results in a dynamical scaling exponent z — 2. This situation must be contrasted to that of a relativistic quantum 
field theory where at tree level (mean field) frequencies and momenta enter with the same power in the propagator 
leading to a dynamical scaling exponent z — 1. Furthermore, critical slowing down is automatically built in the 
phenomenological description, even at tree level as a consequence of dissipative equations of motion^. Clearly this 
is not the case in a relativistic quantum field theory. For a detailed discussion of the differences of non-equilibrium 
dynamical aspects between the time dependent Landau Ginzburg approach and quantum field theory see ref. [^. 

There are important non-equilibrium consequences of slow dynamics near critical points. If the cooling time scale 
is much shorter than the relaxation time scale of long- wavelength fluctuations, these freeze out and undergo spinodal 
instabilities when the temperature falls below the critical^, ^ during continuous (no metastability) phase transitions. 
These instabilities result in the formation of correlated domains that grow in time|^, |^ with a law that in general 
depends on the cooling rate[^. 

In ultrarelativistic heavy ion collisions during the expansion of the quark gluon plasma the critical point for the 
chiral phase transition may be reached. If long wavelength fluctuations freeze out shortly before the transition, the 
ensuing instabilities may lead to distinct event by event observables in the pion distribution as well as in the photon 
spectrum at low energies 

Thus an important aspect of the chiral phase transition is to establish the relaxation time scales of long-wavelength 
fluctuations, and whether critical slowing down and freeze out of long-wavelength fluctuations can ensue. 

In the strict chiral limit with massless up and down quarks, QCD has an SU{2)b, ® SU{2)l symmetry which 
is spontaneously broken to S'[/(2)fl_|_L at the chiral phase transition, the three pions being the Goldstone bosons 
associated with the broken symmetry. It has been argued that the low energy theory that describes the chiral 
phase transition is in the same universality class as the Heisenberg ferromagnet, i.e, the 0(4) linear Sigma model[|j. 
This argument has been used|^ to provide an assessment of the dynamical aspects of low energy QCD based on 
the phenomenological time-dependent Landau Ginzburg approach to dynamical critical phenomena in condensed 
matter]^. While the universality arguments are appealing, a more microscopic understanding of dynamical critical 
phenomena in quantum field theory is needed and has begun to emerge only recently pT| ]. 

In reference a Wilsonian renormalization group extended to flnite temperature was implemented in a scalar 
quartic field theory. In this approach only one loop diagrams enter in the computation of the beta functions, and the 
imaginary part of the self energy, which arises first at two loops order for T > Tc is accounted for by an imaginary 
part in the effective quartic coupling There it is found that the relaxation rate of zero momentum ffuctuations 7 
reveals critical slowing down in the form 7 ~ jT — Td'^ In \T — Tc\ with 1/ ^ 0.53 being the critical exponent for the 
correlation length pO|. 

In reference pil the width of quasiparticles near the critical point has been studied via the large N approximation. 
This study revealed that at high temperature the effective coupling is driven to a (Wilson-Fischer) fixed point, a result 
that is in agreement with the numerical evidence presented in reference While the results in leading order in the 
large N limit found in [pi] hinted at critical slowing down, albeit in a manner different from the numerical evidence 
of reference JTo|, they also hinted at the breakdown of the quasiparticle picture. A conclusion in||ll| is that while the 
large N limit provides a partial resummation of the perturbative expansion, further resummation is needed to fully 
address the relaxation of quasiparticles. 

The large N limit in static critical phenomena presents a similar situation: while it sums the series of bubbles 
replacing the bare vertex by the effective coupling that is driven to the fixed point in the infrared, the self-energy still 
features infrared logarithms that require further resummation |^. Such a resummation is provided by the renormal- 
ization group 

While our motivation for studying dynamical critical phenomena near critical points is driven by the experimental 
program in ultrarelativistic heavy ion collisions to study the QCD phase transitions, the underlying questions are 
more overarching and of a truly interdisciplinary nature. In particular we mention an impressive body of work on 
aspects of quantum phase transitions in condensed matter systems[^ that addresses very similar questions. The work 
in reference|T3] focuses on understanding the static, dynamical and transport properties of low dimensional systems 
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in the quantum regime, in which the frequency and momentum of excitations is uj ; T. 

Our study in this article is complementary to that program in that we focus on the dynamical aspects of long- 
wavelength quasiparticles with lu; p <ti T. As discussed in||l^, |l^ and in detail below this is closer to the classical 
regime. 

The goals: in this article we study the dynamical aspects of quasiparticles near the critical point in a scalar quar- 
tic field theory by implementing a renormalization group program at high temperature. While the renormalization 
group has been generalized to finite temperature in various formulations [l5| , mainly to study critical phenom- 
ena associated with finite temperature phase transitions in field theory, only static aspects were studied with these 
approaches. 

Instead we focus on dynamical aspects, in particular the dispersion relations, and relaxation rates of long- wavelength 
excitations at and near the critical point. Already at the technical level one can see the differences: to understand 
dynamical aspects, in particular relaxation, a consistent treatment of absorptive parts of self-energies is required. 
This aspect is notoriously difficult to implement in a Wilsonian approach in Euclidean field theoryQ. Reference |ic|] 
proposes a method to circumvent this problem but a complete treatment that manifestly includes absorptive parts 
of self-energy contributions is still lacking in this approach. Other approaches using the Euclidean version of the 
renormalization group adapted to finite temperature field theory were restricted to static quantities ||l5|| and in fact, 
as it will be seen in detail below, miss important phenomena that will be at the heart of the results presented here. 

Brief summary of results: Long- wavelength phenomena at high temperature T implies a dimensional reduction 
from the decoupling of Matsubara modes with non-zero frequency ||l^, The coupling in the dimensionally reduced 
theory is AT, where A is the quartic coupling. By dimensional reasons, the perturbative expansion in four space-time 
dimensions is in terms of the dimensionless ratio XT/ fi with the typical momentum scale, which is strongly relevant 
in the infrared. As a result, a perturbative approach to studying long-wavelength phenomena breaks down. This is 
manifest in the breakdown of the quasiparticle picture in naive perturbation theory (see||ll[ and below). 

In 5— e space-time dimensions, the effective coupling in the high temperature, long wavelength limit is g{fJ-) = XTfi^'^. 
We implement an e expansion around five space-time dimensions and a renormalization group resummation program 
at high temperature with T ^ s, ^ near the critical point, with s, fi the typical frequency and momentum scales. We 
analyze the high temperature behavior of the relevant graphs and find that it is dominated for e > by the zero 
Matsubara mode while the sum of the nonzero modes gives subdominant contributions. The effective renormalized 
coupling is driven to an infrared stable fixed point g* = C(e), which for small e allows a consistent perturbative 
expansion near the fixed point. 

An important feature that emerges clearly in this approach, and that has been missed in most other treatments of 
renormalization group at finite temperature, is the anisotropic scaling between spatial and time directions, which is 
manifest in a non-trivial renormalization of the speed of light. This is a consequence of the fact that in the Euclidean 
formulation at finite temperature time is compatified to < r < l/T, thus space and time or momentum and 
frequency play different roles. This results in a novel dynamical critical exponent, z which determines the anisotropic 
scaling. The renormalization group leads to scaling in the infrared region in terms of anomalous dimensions which 
can be computed systematically in the e expansion. In particular, to lowest order in e we find for a scalar theory with 
discrete symmetry z = 1 + + (^(e^), which describes a new universality class for dynamical critical phenomena in 
quantum field theory. All dynamical aspects, such as the relaxation rates, and dispersion relations depend on this 
critical exponent, while the static aspects are completely described by the usual critical exponents. 

We provide a renormalization group analysis of the quasiparticle properties near the critical point, such as their 
dispersion relation and width, complemented with an explicit evaluation to lowest order in the e expansion. The 
main results obtained in this article are the following: at T = Tc we find that the dispersion relation and width 
of quasiparticles of momentum p is Up ^ p^ and Tp ~ {z — l)ujp respectively with a vanishing group velocity of 
quasiparticles in the long wavelength limit highlighting the collective nature of the quasiparticle excitations. For 

T >Tc but \T -Tc\<^T we find uip ~ [l -I- (p^)^^] ^ and Tp ~ (z - 1) ujp j^qj^ with ^ the finite temperature 
correlation length ^ cx |T — Td In the case of 0{N) symmetry we find to lowest order in the epsilon expansion 
that the dynamical exponent z — 1 + e (^^jj^^yi + 0{e'^). 

Critical slowing down emerges near the critical point and in the e expansion Tp/ujp <^ 1 confirming the quasiparticle 
picture. 

We discuss some relevant cases of threshold singularities in which the usual (Breit-Wigner) parametrization of the 
quasiparticle propagator is not available since the real part of the inverse Green's function vanishes at the quasiparticle 
frequency with an anomalous power law. 

In section II we introduce the model and discuss the breakdown of naive perturbation theory. In section III we 
introduce the e expansion and analyze the static case. In section IV the renormalization aspects and the anisotropic 
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scaling is analyzed in detail. Section V presents the renormalization group in the effective, dimensionally reduced 
theory both at and near the critical point. This section contains the bulk of our results which are summarized in section 
VI. Our conclusions and a discussion of potential implications are presented in section VII. The high temperature 
behavior of the relevant diagrams is computed in the appendices. 



II. THE THEORY AND THE NECESSITY FOR RESUMMATION 



The low energy sector of QCD with two massless (up and down) quarks is conjectured to be in the same universality 
class as the 0(4) Heisenberg f erromagnet described by the 0(4) linear sigma model. Furthermore, since we are 
interested in describing the dynamical aspects associated with critical slowing down and freeze out of long-wavelength 
fluctuations just before the chiral phase transition, we focus on T ^ . 

While our motivation for studying critical slowing down stems from the experimental program in ultrarelativistic 
heavy ion collisions, the questions are of a fundamental nature. 

To understand the dynamical aspects near the critical point, we focus on the simpler case of a single scalar field 
theory, and we will recover the case of 0{N) symmetry at the end of the discussion. We thus focus on the theory 
described by the Lagrangian density 

C = l{d,<^>)' + lml<S>'-^<S>' (II.l) 

where the subscripts in the mass and coupling refer to bare quantities. The case of N components in the unbroken phase 
T >Tc differs from the single scalar field by combinatoric factors that change the critical exponents quantitatively. 
These factors will be at the end of the calculation to obtain an estimate of the critical exponents for the 0{N) theory 
in section 



VI 



We are interested in obtaining the relaxation properties of long wavelength excitations near the critical temperature, 
which in this scalar theory Tc ^ \'nij^\/ \/\^ with to_r ; A^, being the renormalizcd mass and coupling. Thus the regime 
of interest for this work is p, lo T ^ Tc with p, oj being the momentum and frequency of the long wavelength 
excitation. As it will become clear below it is convenient to work in the Matsubara representation of finite temperature 
field theory, which is more amenable to the implementation of the (Euclidean) renormalization group. 

In the Matsubara formulation Euclidean time r is compactified in the interval 0<r</3 = l/T whereas space is 
infinite, bosonic fields are periodic in Euclidean time and can be expanded as[E8l Eol pOl 



$(f,T) = _= V / _^0(p',c.„)e-*""^+^P"-" (n.2) 




(II.3) 

Thus we see that while the spatial momentum is a continuum variable, the Matsubara frequencies are discrete as a 
consequence of the compactification of Euclidean time. This feature of Euclidean field theory at finite temperature 
will be seen to lead to anisotropic rescaling between space and time and therefore, as it will be clear below, new 
dynamical critical exponents. Anticipating anisotropic rescaling, we then introduce the bare speed of propagation wo 
of excitations in the medium by writing the Euclidean Lagrangian in the form 

= \^^^ + ^(^*)' + ^^^'(^) + ^ + $2 + ^ $4 (11.4) 

where we have introduced the effective renormalized temperature dependent mass M(T) and the counterterms, in 
particular the mass counterterm Srn^iT) = —M'^(T) — tuq is adjusted order by order in perturbation theory so that 
the inverse two point function obeys 

r(2)(p = 0;c^„ = o) = M2(r) (n.5) 

The critical point is defined with (the inverse susceptibility) M'^iT) = 0. We will begin our study by focusing our 
attention on the critical theory for which M{T) — 0. We will later consider the theory near the critical point but in 
a regime in which M{T) <C T ~ Tc- Thus, the general regime to be studied is p, M{T) <^T ^ Tc- 



A. Infrared behavior of the critical theory: static limit in three space dimensions 



In order to highlight the nature of the infrared behavior when p,uj<^T we focus first on the critical theory in the 
static limit, when the Matsubara frequencies of all external legs in the n-point functions vanish. For the purpose of 
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understanding the nature of the infrared physics in the static hmit we will set iiq = 1 in this and next section, and we 



will recover this variable, the speed of light, when we study the dynamics in section IV 



1. The scattering amplitude in D — 3 

We consider first the 2 — > 2 scattering amplitude, or four point function, to one-loop order in three spatial dimen- 
sions. The full expression is given by 

r''*'(pi, Si,P2,S2,P3, S3,P4,S4) = - Xq + [H {pi + P2, + S2) + H {pi + P3, Si + S3) 

+ H{pi+p4,si+Si)]+0{\l) , (II.6) 

where = 27r T , 1 < « < 4 and to,; G Z, 

" 2" 2 / (27r)3 [g2 + {2tt T n)^] [{q + p)^ + (2^ {n + mf] ' ^^^''^^ 

s = 2'kT m. Since the external momentum p <C T it is clear from the above expression that the dominant infrared 
behavior of H{j)^ 0) is determined by the zero Matsubara frequency in the sum. As it will be explicitly shown below 
the contribution from the non-zero Matsubara frequencies will introduce a renormalization of the bare coupling which 
in the limit T ^ p\s independent of the external momentum (this will be seen explicitly in the next section). Keeping 
only the zero internal Matsubara frequency and carrying out the three dimensional integral explicity we find 

i^^r(p,0)=-^ (II.8) 

Ibp 

Thus, defining the effective coupling constant at the symmetric point pi = Pi where 

A • P, = (4 5,, - 1) ^ (II.9) 
in the static limit one finds that in the infrared limit ^i/T <^ 1 



'^effi^') — Aq 



o{xi) (11.10) 



Two important features transpire from this expression: i) the factor T//z can be explained by dimensional arguments: 
in the Matsubara formulation for each loop there is a factor T from the sum over internal Matsubara frequencies. The 
infrared behavior for /z <C T is obtained by considering only the zero internal Matsubara frequency in the loop. This 
integral has only one scale and since A is dimensionless in three spatial dimensions the one-loop contribution must be 
proportional to T/ ^. A similar argument shows that for a diagram with m internal loops, and transferred momentum 
scale /i there will be a power T™ from the Matsubara sums, the infrared behavior is obtained by the contribution 
with all the internal Matsubara frequencies equal to zero, which by dimensional power counting must be of the form 
(T/zi)™. Therefore a diagram with m-internal loops will contribute to the scattering amplitude with A(AT//i)™. In 
taking only the zero internal Matsubara frequency we are assuming that the internal loop momenta arc cutoff at a 
scale < T. 

Thus, at the critical point the most important infrared behavior is that of the dimensionally reduced three di- 
mensional theory [|T6| [T^ . The reason for this dimensional reduction is clear: at finite temperature T the Euclidean 
time is compactified to a cylinder of radius L = 1/T for transferred momenta /i the spatial resolution is on distances 
d ^ Therefore ioi fi <^ T ^ d ^ L thus the compactification radius is effectively zero insofar as the long distance 
(infrared) physics is concerned. 

We will study below the contribution from the non-zero Matsubara frequencies, ii) for a transferred momentum scale 
/i perturbation theory breaks down for /i <C AT since the contribution from higher orders is of the form A(Ar//x)™. 
This suggests that a resummation scheme is needed to study the infrared limit. This situation is similar to that in 
critical phenomena where infrared divergences must be summed and the renormalization group provides a consistent 
and systematic resummation procedure. We can obtain a hint of how to implement the renormalization group in finite 
temperature field theory in the limit when T is much larger than any other scale (masses, momenta and frequency) by 
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realizing that, from the argument pres ented above, the perturbative expansion is actually in terms of the dimensionless 
coupling gQ — XT /pi. Therefore from ( ILlOj ) we can write 



9effifJ^) = 90 



, 3 



Oigl) 



We can improve the scattering amplitude via the renormalization group by considering the RG-/3 function 

3 



-9eff + ^9iff + 0{g, 



(11.11) 



(11.12) 



The first term (with the minus sign) just displays the scaling dimension (for fixed AT) of the effective coupling, that 
this dimension is —1 is a consequence of the dimensional reduction since AT is the effective dimensionful coupling of 
the three dimensional theory. 

We thus see that the renormalization group improved coupling runs to the infrared fixed point 



9 



(11.13) 



as the momentum scale — > 0. Comparing with the renormalization group beta function of critical phenomena[^ |2^, 
p^ , p3| , p^ , we see that this is the Wilson-Fischer fixed point in three dimensions, again revealing the dimensional 
reduction of the low energy theory. The resummation of the effective coupling and the fixed point structure can 
also be understood in the large N limit pl]|. As described in|l^ the large N limit can be obtained by replacing the 
interaction in the Lagrangian density forpGl 



2N\ I 



(11.14) 



with $ = ((/)i , • ■ ■ (^jv), and the form of the quartic coupling has been chosen for consistency with the notation 
of reference [|ll|. The leading order in the large TV limit for the scattering amplitude is obtained by summing the 
geometric series of one-loop bubbles in the s-channel (only this channel out of the three contribute to leading order 
in the large N limit), each one proportional to N which is the number of fields in the loop. As a result one finds that 
the effective scattering amplitude at a momentum transfer is given bypl| 



A 



AT 
4^ 



Thus, introducing the dimensionless effective coupling gf,ff{pi) — Xef f{fJ.)(T/ pi) one finds that 

Iim5e//(M) = 4 



(11.15) 



(11.16) 



i.e, the effective coupling constant goes to the three dimensional fixed point |nT| 



2. The two point function in D = 3: 



The two point function in the static limit is given by 

r(2)(p,0) = / + 5m2(T) 



E(p,0)-KO(A3) 



(11.17) 



where E(p, 0) stands for the two loop sunset diagram at zero external Matsubara frequency and the counterterm 
5vn?[T) will cancel the momentum independent, but temperature dependent parts of the self-energy. The two loop 
self energy for external momentum p and Matsubara frequency LUm — iirTm is given by 



E(p, w„ 



A2T2 



E 

LjeZ 



<pk 



(2^)3 (2^)3 



{p + k + q)'^ + {UJI + LUj + LOmY 



(11.18) 



with LOj — 27: jT. The static limit is obtained by setting LOm = (m = 0). In this limit the dominant contribution in 
the infrared for T ^ p arises from the term I — j — in the sum. The terms 1^0; j 7^ for which we can take 
p = (since p <^T) will be cancelled by the counterterm. A straightforward calculation leads to 



T(^\pMm=P^ 



1 



In 



12(47r)V Vm' 



P 



(11.19) 
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where /i is a renormalization scale. This expression clearly reveals the effective coupling XT /p which becomes very 
large in the limit p <^ XT. Clearly we need to implement a resummation scheme that will effectively replace the 
bare dimensionless coupling constant by the effective coupling that goes to a fixed point in the long- wavelength limit, 
and also ensure that at this fixed point the effective coupling is small so that perturbation theory near this fixed 
point is reliable. This is precisely what the renormalization group combined with the e-expansion achieves in critical 
phenomenal, H, H, H, H. 



B. Dynamics in D = 3 

The two- loop contribution to the self-energy for w 7^ is obtained from a dispersive representation of the self-energy 
in terms of the spectral density 



■K V — LO — iO+ 



(11.20) 



The spectral density p(p, v) has been obtained in reference pH in the high temperature limit pTf in D = 3. Using the 
expression given in |ll[| for the spectral density at two loops in the high temperature limit, and after some lengthy but 
straightforward algebra, we find 



12 V 47r 



Q{\u\~p) + l^Q[p-\v\) 
P 



(11.21) 



Carrying out the dispersive integral (11.20) and subtracting off the terms that are independent of p, uj that are absorbed 
by the counterterm we find 



' 12 V 47r 



2 r 



In 



p2 - (cj + iO+f 

7 



LU /' UJ + iO^ — p 
— In 
P 



uj + iO+ +p 



(11.22) 



Clearly, the static limit u; of the self-energy coincides with ([1.19). The two point function is therefore given by 

1 /XT' 



r(2)(p. 



2 2 

-- P - bJ + 



12 V 47r 
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f 


In 






LO — p 








^2 


P 


UJ +p 





ip{p,U}) 



(11.23) 



with p(p,uj) given by ( [1.21 ). 

There are several features of this expression that are noteworthy: 

• It is clear that for AT ^ p;ct; the two loop contribution is much larger than the tree level term p^ — this 
already signals the breakdown of perturbation theory in the high temperature regime when AT ^ p\uj in the 
dynamical case. 



Consider the real part of the two point function as cj — > p i.e, near the mass shell. 



ReT^^\p,uj^p)^2p^ (1-^ 



\^-p\ 
P 



1 (^\\^C^ 

12 \'i-Kp) V /i 



(11.24) 



This expression reveals that uj — p \& not the position of the mass shell of the (quasi)-particle. The coefficient 
of (1 — uj/p) hints at wave function renormalization but the fact that the two point function does not vanish at 
this point prevents such identification. Furthermore, we see that the term that does not vanish at ui — p hints 
at a momentum dependent shift of the position of the pole, i.e, a correction to the dispersion relation. However 
for XT /p 3> 1 both contributions are non-perturbatively large and the analysis is untrustworthy. 

• Now consider the width of the quasiparticle. 



7p 



ImY,[p,uj = p) T^P ( XT 



2p 



12 \4ttp 



(11.25) 



so that 7p/p ^ 1 for XT/p ^ 1. This signals the breakdown of the quasiparticle picture. 
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A similar analysis reveals the breakdown of perturbation theory away, but near the critical point with \T — Tc\ <C 
T . The imaginary part of the two-loops self-energy aX p — Q and in terms of the temperature dependent mass 
rrifiiT) can be obtained straighforwardly in three spatial dimensions'^ in the limit T 3> m{T). It is found 
to be|ll|, |l|] 

ImE(2)(p'^ 0,w = TOfl(T)) cx A^T^ (11.26) 

consequently, at two-loops order the width of the zero momentum quasiparticlc in three spatial dimensions is 
given by|l|, |l| 

r (X » mfl(r) (11.27) 



• This behavior is different from that of gauge theories at high temperature where low order fermion or gauge 
boson loops are infrared safe and determined by the hard thermal loop contributions . This is so because for 
fermions there is no zero Matsubara frequency while in the case of gauge bosons, the vertices are momentum 
dependent. While the second term inside the bracket in expression ( [L23| ) is determined by Landau damping is 
infrared finite and is similar to the leading contribution in the hard thermal loop program jl^, the first term 
arises from the three particle cut. The dependence of this term on the renormalization scale /i arises from the 
subtraction of the mass term at the critical point and reveals the infrared behavior. Furthermore, in the hard 
thermal loop program |l9| ^ one finds thermal masses of order gT with g the gauge coupling, and widths of 
order g^T (up to lo garithms) so that Tp/uip ^ 1 in the weak coupling limit, while in the scalar theory under 
consideration (11.27) suggests that Fp/wp ^ 1 in naive perturbation theory. 



• We note at thi s stag e that the high temperature limit p, w ^ T of the self-energy calculated from the spectral 
r epres entation ( 11.20 ) can be directly obtained by computing I](p, a;„i) in the Matsubara representation given by 
([L18) by setting the internal Matsubara frequencies loi = loj = and analytically continuing cOm ~* ~i^ + O"*". 

We highlight this observation since it will be the basis of further analysis in what follows: the high temperature 
limit of the self energy p,Lu ^ T can be obtained by setting the internal Matsubara frequencies to zero and 
analytically continuing in the external Matsubara frequency, i.e. 



S(p,w) 



d^q d^k 



6 J (2^)3 (27r)3^2p (^^+k + q)^+u}l 



(n.28) 



We provide one and two loops examples of this statement in appendices ^ and ^ and formal proof of this 
statement to one loop order in appendix |c[ 

The result for the width of the quasiparticlc at two loops was anticipated in references [0 This width is 
purely classical since the product AT is independent of hfVi^ . This result for the damping rate of long wavelength 
quasiparticles in the critical theory is in striking contrast with that for T » Tc which has been studied in detail 
inp9|, |30[ |. For T » the thermal mass is ruth oc \/ATp9| while the two loop contribution to the imaginary 
part of the self energy for T » Tc,p is still oc A^T^. Thus, the damping rate of long wavelength excitations 
is 7 oc A^/^T <^ A^/^T in the weak coupling limit. Therefore for T » Tc long wavelength excitations are true 
weakly coupled quasiparticles with narrow widths. 



III. THE e-EXPANSION: STATIC CASE 



The analysis of the previous section points out that naive perturbation theory at finite temperature breaks down 
at the critical point for momenta <C AT. The reason for this breakdown, as revealed by the analysis of the previous 
section, is the following: in four space-time dimensions the quartic coupling is dimensionless, however each loop 
diagram in the perturbative expansion has a factor T from the sum over the Matsubara frequencies. After performing 
the renormalization of the mass including the finite temperature corrections and setting the theory at (or near) the 
critical point, the effective expansion parameter for long- wavelength correlation functions is AT which has dimensions 
of momentum. If a given diagram has a momentum transfer scale /x the effective dimensionless expansion parameter 
is therefore XT/ fi which becomes very large for fi <^ AT, i.e, the effective coupling is strongly relevant in the infrared. 
The analysis based on the RG beta function ([1.12, [1.13) suggests that the effective coupling g — XT/ ji is driven to the 
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three dimensional (Wilson-Fischer) fixed point in the infrared, obviously a consequence of the dimension al reductio n 
in the high temperature limit. This is confirmed by the large N resummation of the scattering amplitude (11.15 11.16). 
// the value of the coupling at the fixed point is <C 1 then a perturbative expansion near the fixed point would be 
reliable, however the value of the coupling at the fixed point is g* ^ 0{1) which of course is a consequence of the fact 
that for fixed AT the effective coupling scales with dimension of inverse momentum in the infrared. This situation 
is the same as in critical phenomena for theories that are superrenormalizable, in which the infrared divergences are 
severe. 

The remedy in critical phenomena is to study the perturbative series via the e-expansion wherein the value of the 



coupling at the fixed point is C(e)y, 21, |2^, ^ and sum the perturbative series via the renormalization group. 
We now implement this program in the high temperature limit. 

In five space-time dimensions the quartic coupling A has canonical dimension of inverse momentum, therefore the 
product AT that occurs in the perturbative expansion in the dimensionally reduced low energy theory is dimensionless. 
Then in a perturbative expansion at (or very near) the critical point we expect that infrared divergences will be 
manifest in the form of logarithms of the momentum scale in the loop. This implies that the effective coupling is 
marginal. Considering the theory in 4 — e spatial dimensions and one Euclidean (compactified) time dimension the 
effective coupling of the dimensionally reduced theory, AT, has dimensions of with /i being a momentum scale. 
Therefore the effective dimensionless coupling for diagrams with a transferred momentum scale /i is g{^) — AT/i^*^. 
Thus, for fixed T the scaling dimension of this effective coupling is — e, hence we expect a non-trivial fixed point at 
which the coupling g* ^ 0{e). 

Therefore for e ^ 1 we can perform a systematic perturbative expansion near the fixed point. This is the spirit of 
the e expansion in critical phenomena which, when combined with a resummation of the perturbative series via the 
renormalizationgroup has provided a spectacular quantitative and qualitative understanding of critical phenomena|^. 

While dimensional regularization and the e expansion have been used to study the dimensionally reduced high 
temperature theory insofar as thermodynamic quantities is concerned, i.e, static phenomena ^ we emphasize 
that our focus is to study dynamics at and near the critical point, which is fundamentally different from the studies 
of static phenomena in these references. 

As a prelude to the study of the dynamics, we now revisit the scattering amplitude at one loop level and the self- 
energy at two loops level in 4 — e spatial dimensions at high temperature in the static limit. The one loop self-energy 
is momentum independent and is absorbed in the definition of the thermal mass [p9| which is set to zero at the critical 
point. The are two main purposes of this exercise: the first is to quantify the role of the higher Matsubara modes 
and secondly to obtain a guide for the infrared running of the coupling constant. 



A. Scattering amplitude: 



The one loop contribution in the static limit, i.e, when the external Matsubara frequencies are zero is given by 



H{p,Q) 



T 



E 



1 



[AttY [g2 + (27r T nf] [{q + p)^ + {2n T 



(III.1) 



In the high temperature limit the nonzero Matsubara terms give subdominant contributions. This property can be 
argued in different manners. For example, the / ^ terms in eq.( III.l ) can be interpreted as Feynman diagrams in 
d — A — e dimensions with mass {2n T l)^. Such contributions are negligible for p ^ T[^, 
We find for high temperatures (see Appendix ^ ), 



F(p,0)^>^-i/..(p,0)-T-^^(^ + ^) ^^'-"'^P' 



192 TT^+f 



rp2 



l + O 



El 

rp2 



(III.2) 



where Hasiip, ) stands for the I — contribution to the sum (ill.l ) plus the dominant high temperature limit of the 
sum over the / terms, which is obtained by setting p — 0. Separating the I = mode and setting p = in the 
contribution from the sum over I ^ we find 



Has^{p, 0) = 



Tp- 



2(47r) 



r^(i-f) 
r(2-e) 



with ^ Riemann's zeta function which has the following properties! 

1 1 

2 



C(o) 



lim C(e) 



4^ 

p2 



+ 1 



C(e) 



(III.3) 



(III.4) 



10 



where 7 = 0.577216 . . . stands for the Euler-Mascheroni constant. 
Near four space dimensions e we find 



Has^{p, 0) 



T 



2(47r) 



log 



p 



4T2 



0(e) 



(III.5) 



There is no pole in e and the argument of the logarithm reveals that T acts as an ultraviolet cutoff. The reason 
that there is no pole in epsilon as e — > is that the poles in e should be independent of temperature and should be 
those of the zero temperature theory. However, in dimensional regularization one loop integrals have no poles in odd 
space-time dimensions [p4| . On the other hand near three space dimensions e — > 1 we find 



Hasr{p,0) 



T 

IQp 



1 



— ^In- 



1 



(47r)2 (e- 1) 



(III.6) 



where the pole term at e = 1 corresponds to the usual coupling constant renormalization. This divergent term is 
temperature-independent, as expected, the ln(r) is reminiscent of an u pper momentum cutoff for the high temperature 
limit. The first term is precisely what we obtained in equation (II. 8) from setting the internal Matsubara frequency 



to zero, i.e, the result of the dimensionally reduced theory. After subtracting the pole near three space dimensions, 
the first term gives the leading infrared contribution in the limit T /p 3> 1 whereas the logarithm is subleading. 

This expression coincides with that given in references p6|, p^ . In these references the four dimensional high 
temperature static theory was studied and a systematic analysis of Feynman diagrams in the dimensionally reduced 
theory (three dimensions) was performed. The e — 1 in our expressions should be mapped onto the 2e for comparison 
with the results in these references. 

For e > we can neglect the terms of the form {T /p)~^ in eqn. ( 111.3 ) in the limit T /p ^ 1. And for 1 ^ e > we 
find that the static scattering amplitude at a symmetric point ( p^I.9| ) in the limit in which the temperature is much 
larger than the external momentum scales is given by 



r(4)(p, = F„o) = -Ae//(/z,r) = -A 



3 \Tii- 



(47r) 



(III.7) 



the factor T arising from the Matsubara sum is such that AT has dimensions of (momemtum)*^ so that in 5 — e space 
time dimensions AT/i"*^ is dimensionless. Thus, introducing the dimensionless renornialized coupling 



(47r)i 



we find 



-£5+35^ + 0(3^) 



Therefore this effective coupling in the infrared limit is driven to a non-trivial fixed point 

5=0 



(III.^ 



(III.9) 



(III.IO) 



Hence for e <C 1 the fixed point theory can be studied perturbatively. This of course is the basis of the e expansion 
in critical phenomenal^, |2^, |^, |2^ and will be the important point upon which our analysis will hinge. 



B. Two loops self-energy: 



As mentioned above the one loop contribution to the self energy is momentum independent and absorbed into the 
definition of the thermal mass. The two loop contribution in the static limit in d spatial dimensions is given by 



E(2)(fc,0) 



A2T2 



EE 

/ rii 



d<^p d'^q 



1 



{2ttY (27r)'^ + c^2][p2 + ^2J[(p + 5 + fc)2 + (^^ + 



(IH.U) 
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We now introduce two Feynman parameters, separate the Z = m = term from the Matsubara sums and take T ^ p 
in the sums with Z, m ^ to find 



rpw ^ 9^r(-l + e)r3(l-f) „ 



6 



r(3-f) 





— e 

















C{d) 



(III.12) 



with C{d) only depending on the dimensionahty. The second term oc T** ^"^ does not depend on the momentum and 
is therefore cancelled by the mass counterterm which defines the critical theory. Therefore for e > but small we find 



r(=^)(fc,0) _ [S(A:,0) -E(0,0)] 



'+12^-12 



We introduce the wave function renormalization in dimensional regularization by the usual relation 
and choose 

Therefore the renormalized two point function in the static, high temperature limit is given by 



l-Aln 
12 



(III. 13) 



(in. 14) 



(III.15) 



(III.16) 



The infrared behavior is obtained by resumming the perturbative series via the renormalization group |2^, 
p4, Ea| which leads to the scaling form of the two point function in the infrared limit 



r^^'(fc) cx fc^-"? 



(III.17) 



with 



6 54 



(III. 18) 



This is the anomalous dimension to lowest order in e 14, G^, \22, 23 



H P5 



C. The strategy: 

The analysis of the static case above has highlighted several important features of the infrared behavior near the 
critical point, which determines the strategy for studying the dynamical case: a:) the infrared behavior in the limit 
when T 3> p with p the typical momentum of the Feynman diagram is determined by the dimensionally reduced 
theory obtained by setting the internal Matsubara frequencies to zero, b:) naive perturbation theory breaks down 
in three space dimensions because the dimensionless coupling is XT/ fi with fj, the external momentum scale in the 
Feynman diagram, while a large N or renormalization group resummation suggests a non-trivial infrared stable fixed 
point, the coupling at this fixed point is of C(l). c:) Just as in critical phenomena the perturbative expansion can be 
systematically controlled in an e expansion around four spatial dimensions corresponding to a theory dimensionally 
reduced from 5 — e space time dimensions. The effective dimensionless coupling of the dimensionally reduced theory 
(four dimensional) is AT which is marginal. This combination is independent of ?i-this can be seen by restoring powers 
of /i A — > hX; T — !■ T /h - so that the low energy, dimensionally reduced theory is classical. In 4 — e spatial dimensions 
the effective dimensionless coupling g — XTp,^'^ is driven to the infrared stable Wilson-Fischer fixed point of 0{e) by 
the renormalization group trajectories. Thus, the strategy to follow becomes clear: we will now study the dynamics 
by including the contribution from the external Matsubara frequency, focusing on the infrared behavior for p,uj <ti T 
near the critical point in a systematic e expansion around four spatial dimensions. 

We note that the theory in 5 — e dimension is formally non-renormalizable in the ultraviolet, however this is 
irrelevant for the infrared which is the region of interest here. The analysis provided above in the limit e — > clearly 
shows that near five space-time dimensions there are no poles in dimensional regularization in one loop diagrams as 
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expected [plf. The potential poles are replaced by ln(T). The low-energy theory must be understood with a cutoff of 
0(T) and the dimcnsionally regularized integrals in 5 space time dimensions clearly display this cutoff in the argument 
of logarithms. The long wavelength fi/T <C 1 an the e limits do not commute: keeping the subleading terms 
in the high temperature limit and taking e —^ results in that poles in e actually translate into logarithms of the 
cutoff T, on the other hand, keeping e > and small the T/fi — > oo limit can be taken and the subleading high 
temperature corrections vanish. Clearly it is the latter limit the one that has physical relevance, since eventually 
we are interested in studying the infrared behavior of the physical theory in three space dimensions. Hence in what 
follows we consider the long-wavelength limit for e > but small and approach the physical dimensionality e — > 1 in 
a consistent e expansion improved via the renormalization group. This is the strategy in classical critical phenomena 
as well where for e > and small the ultraviolet cutoff can be taken to infinity. 

At this stage it is important to highlight the difference between the main focus of this work and that in references |l6| 
[iTt . The work of references ||l^, |l^ studies the static limit of the dimcnsionally reduced theory near three spatial 
dimensions arising from the high temperature limit of a four dimensional Euclidean theory compactified in the time 
direction. In contrast we here focus on studying the dynamics in the limit when p,Lu <^ T which as emphasized by 
the analysis above will be studied in an e expansion in a dimcnsionally reduced theory near four space dimensions. 

The limit of phy sical interest e — s- 1 must be studied by improving the perturbative expansion via the renormalization 
group [|, m, 11 H, m H and eventually by other non-perturbative resummation methods, such as Pade approximants 
or Borel resummation that will extend the regime of validity of the e expansion p^. 



IV. DYNAMICS NEAR THE CRITICAL POINT: 



We now turn to the dynamics. Our main goal is to study the feasibility of a quasiparticle description of low energy 
excitations at and near the critical point. Of particular interest is the dispersion relation as well as the damping 
rates of these excitations. This information is contained in the two point function r(^^(p, u;„) which is the inverse 
propagator, analytically continued to Un —iw + 0+. The region of interest is p,uj ^ T and if the theory is (slightly) 
away from the critical point M{T) ^ T as well. In principle for a fixed (jJ„ or fixed external Matsubara frequencies 
in the external legs of n — point functions, one must perform the sum over the internal Matsubara first and then 
take the analytic continuation. However, as was shown above in detail in the static case, the most infrared singular 
contribution arises from setting the internal Matsubara frequencies to zero. That this is also the case in the dynamics 
can be seen by considering a diagram with m— internal lines, rerouting the external Matsubara frequency through 
one of the lines. All of the lines are equivalent since rerouting the external Matsubara frequency corresponds to a shift 
in one of the sums. The other m — 1 contain propagators in which the internal Matsubara frequency acts as a mass of 
0{2ttI T). These are the superheavy modes in the description of references p^, p7| pl| , |3^ . The contribution that is 
dominant in the infrared is from the region of loop momenta <C T which is largest when the mass of the propagator 
is zero, i.e, the zero Matsubara frequency. Keeping non-zero Matsubara frequencies in any of the rn — 1 legs will lead 
to subleading contributions in the limit p, w <C T. 

Once the internal Matsubara frequencies had been set to zero we can analytically continue the external Matsubara 
frequency to a continuous Euclidean variable aj„ — s- s to obtain the Euclidean two point function. The dispersion 
relation and damping rate are obtained by further analytical continuation s —iu + 0^. 



As anticipated in section(II) because Euclidean time is comp actif ied and plays a different role than the spatial 



dimensions, we must consider the anisotropic Lagrangian density (II.4) which includes the velocity of light multiplying 
the derivatives with respect to Euclidean time. If this velocity of light is simply a constant it can be reabsorbed into 
a trivial redefinition of the time variable. However, as it will become clear below, this velocity of light acquires a 
non-trivial renormalization as a consequence of the anisotropy between space and time directions at finite temperature 
and will run with the renormalization group. Thus, the Euclidean propagator is generalized to 

G(fc,w™) = — ^ (IV.l) 
fc H — T- 



A. The scattering amplitude 



We begin by studying the scattering amplitude, now as a function of external momenta and frequencies. The one 
loop contribution is determined by the function H{p,s) given by eqn. (H.7), which for p, s <C T is given by (see 
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Appendix) 



p2 + :^(l-e) 



l + O 



„2 



y2 



(IV.2) 



where Hasi{p, s) stands for the I = contribution to the sum (11.7 ) plus the dominant high temperature hmit of the 
sum over the / ^ terms which is obtained by setting p, s = 



rr , ^ T f d'^q 
gmequationHasi(P,s) = — 



r(i + i) C{e) ^1 



{2n) 



iq + p)^ + 



This integral is computed in the appendix A with the result 



Hast{p, s) 



r(§-i) 

2(47r)^"* 



T\-+p- 



+ 



r(i + i) C(e) ^1 



(IV.3) 



(IV.4) 



where -^(0, 5; c; z) stands for the hypergeometric function. For e > and p, ^ <^T we can neglect the second term, 
since it is proportional to T^~^ <^ T (p^ + . We note that the infrared dominant contribution can be written in 

the form Tp~'^!F{-^^:^) the factor T thus combines with the coupling A to give the effective coupling of dimension /j,*^ 
in d = 4 — e spatial dimensions. 



For e > but small and neglecting the second term in ([V.4) can be expanded in e leading to 



\{)Hasi{p,s) 



9ip) 



- 1 



9 9 

VqP'' 



In 



+ p'^Vq 



In' 



VqP'^ 



In 471 + 2-7 + 0(6) 



where we introduced the dimensionless bare coupling 

5(m) 



(47r)i 



(IV.5) 



(IV.6) 



We remark that one cannot take e — > in this expression since in this limit the pole is actually cancelled by the 
second term in (IV.4) above. As emphasized above, this expression must be understood for e > but small so that 
the contributions of the form {T/s, T/p)~'^ for T ::§> s,p. Therefore, the expression above must be understood in 
the sense that i): T ^ s,p with fixed e > and ii): e ^ 1 and the resulting expressions have a Laurent expansion for 
small e. 



B. The self energy at two loops 

Neglecting the contribution from the non-zero Matsubara frequencies which will be absorbed by the mass countert- 
erm in the definition of the critical temperature (or M^(T) away from the critical point) and also neglecting terms 
that vanish in the limit T 3> p, s the dominant contribution in the infrared to the two loops self-energy is 



E(p, s) 



(Pk 



d'^q 

(2^7 (27r)%2p 



1 



(g+/c+p)2 + ^ 



AgT2 r2(f -i)r(3-d) /•! 



6(47r) 



r(fi- 2) 



n d-3 



xp -\ 2 



(IV.7) 



While this expression can be written in terms of hypergeometric functions, it is more convenient to expand it in e 
with the result that 



r(2)(p. 



9M 
6e 



'"p^ 2 p^ , I p" 
\ ^ — e — In 



V, 







(IV.8) 



with g(/i) given by eqn. (IV.G) 
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Renormalization 



The forms of the two and fo ur p oint functions immediately suggests a renormalization sche me ak in to the familiar 
one used in critical phenomena pl[ |2^ , p3| , [2^ with one important difference: we see from eqn. ( IV.8| ) that the velocity 



of light vq must also be renormalized. The wave function renormalization is introduced as usual via 
the renormalized mass as a function of temperature is defined as 



(IV.9) 



rg^(o,o) = Af2(r) 



(IV.IO) 



this definition, however, defines the inverse susceptibility or correlation length, rather than the pole mass, the critical 
theory is defined by M'^{T) = 0. Coupling constant and velocity of light renormalization are achieved by 



Z\ Ao 



2 2 Zv 





The renormalization conditions that determine the constants Z,p, Zx, 



(IV.ll) 
(IV.12) 



dp'^ 



T'^^\p,^P,;s^O)^-Xi 



(IV. 13) 



Consistently with the e expansion we choose the renormalization constants Z^,Zx,Zv in the minimal subtraction 
scheme to lowest order, since keeping higher powers of the coupling or e results in higher order corrections in the e 
expansion. 

To lowest order, one loop for the 4-point function and two loops for the two point function, we find from the results 



([V.5) and (IV 



^0 



35(y") 



9B. = g{tj) 



9R 



XrT fi 
(47r)5 
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12e 

2 



g 

3e2 



Zy = 1 

Thus, the renormalized two point function reads 



Tf{p,s)^p^ 



hi;'" 




52 
+ — 

















(IV. 14) 
(IV. 15) 

(IV. 16) 



THE RENORMALIZATION GROUP: 



Before we embark on the resummation program via the renormalization group, it is important to highlight two 
important features 
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• The contributions that are dominant in the infrared in the Umit T ^ p, s correspond to the terms with internal 
Matsubara frequencies equal to zero, the non-zero Matsubara frequencies give subleading contributions for e > 0. 
This in turn results in that the dependence on temperature is solely through the effective coupling g = XTji"'^. 
This can be seen from the fact that each loop has a factor T from the Matsubara sum over the internal loop 
frequencies as well as one power of the coupling constant A, by dimensional reasons the dimensionless coupling 
is obtained by multiplying by /i"'^. 

• The velocity of light v always enters in the form s/v since this is the form that enters in the propagators and 
the renormalization conditions above. 



A. The critical point: 



The bare n-point functions are independent of the renormalization scale and this independence leads to the 
renormalization group equations (we now suppress the subscript R understanding that all quantities are renormalized) 



d_ 



d 



with 



To lowest order wc find 



d N 



' [Pi,—;P2, — ■ ■ ■ ,pn; — ■,g,^i]^o 



dv 



Xa,T,vo 



9/1 



f3g 
Pv 



3e 



v + Oig',g'e) 



Oia^gh) 



(V.l) 

(V.2) 
(V.3) 
(V.4) 

(V.5) 
(V.6) 

(V.7) 



While we can write down the general solution of the RG equation ( V.l ) for an arbitrary iV-point function, our focus 
is to understand the quasiparticle structure which is obtained from T^. 
Since F^^^ has dimension two, it follows that 



r(2)(p,^,5,M) = M^<f (-,-,.9 



(V.8) 



therefore 



This scaling property then leads to 



d d ^ 



r(2)(e*p,l^,5,/.) = 

V 



(V.9) 



(V.IO) 



which combined with the RG equation (V.l) leads to the following equation that determines the scaling properties of 
the two point function 



(v.ii) 
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The solution of this equation is standard ||2l|, g2| g3| 



r(2) 



with 



dgit) 



(V.12) 



dt 
dvjt) 
dt 



(V.13) 



lit) = 7(5(i),«W) 

As t ^ — oo i.e, the momentum and frequency are scaled towards the infrared we see from the RG [3 function ( |V.5| ) 
that coupling is driven to its fixed point 



, lim 9{t) =9* = ^ + 0{^) 



which in turn impHes that 



(V.14) 



lim 7(i) = ri= — 
t^-co 54 



lim v{t) = v{0) e^^-^'* 



t — > — oo 

where we introduced the new dynamical critical exponent 

1 



Vt 



2g* 



3e 



Therefore, in the asymptotic infrared limit we find that 



e(i-^)* 



,9 ,M 



It is convenient to re-define pe* = P ; se* = S* to find 



r(2) ^ 



,9,f^ 



6* 



and finally writing P =^ fj,e^ and using the property (V.8) we find the scaling form in the infrared limit 



.(2) 



P, -,5,Ai 

V 



1 2-n 



^ = 



1 — Z IDZ 



(V.15) 
(V.16) 



(V.17) 



(V.18) 



(V.19) 



(V.20) 



The solution of the RG equation clearly shows that the two point function in the infrared limi t is a s caling function of 
the ratio s/v{ii)ii^^^ p^ hig hlighting the role of the new dynamical exponents z given by eqn. (V.17) with r\t to lowest 
order given by eqn. ( V.l?] ). 

The emergence of the new dynamical exponent z is a consequence of the anisotropic renormalization between mo- 
mentum and frequency, or space and time manifest in the renormalization of the speed of light. This novel phenomenon 
can be traced back to the different role played by time (compactified) and space in the Euclidean formulation at finite 
temperature. A similar anisotropic rescaling emerges in a different context: a Heisenberg ferromagnet with correlated 
impurities with similar renormalization group results. 

While the formal solution does not yield the function <I>, we can find it by matching to the lowest order perturbative 
expansion ( pV.16D when the coupling is at th e non- trivial fixed point. From the form of the of the perturbative 
renormalized two point function given by eqn. ( IV.16| ) and assuming the exponentiation of the leading logarithms via 
the renormalization group near the non-trivial fixed point 



V'-I\v,s-9*)^V^ 
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(V.21) 
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which can be immediately written in the scaUng form 



2-ry 1, 

~ 1 + -[rit~V) 

2-rit 2 



(V.22) 



Clearly this form coincides with the scaling solution of the renormalization group and the perturbative expansion in 
the regime in which it is valid. We note, however, that in the computation of rjt we have neglected contributions 
to the renormalization of the velocity v of 0{g^/e) which would appear at next order and would lead to an 0{£^) 
contribution to rjt, which is of the same order as rj in z. Thus consistently we must neglect the contribution of 77 to 
the dynamical exponent z which to lowest order is therefore 



= 1 



(V.23) 



1. Quasiparticles and critical slowing down: 

The quasiparticle structure of the theory is obtained from the Green 's function G-i(p,a;) =T'-'^^p,s = -iuj+0+,g*). 
In particular the dispersion relation and the width of the quasiparticle are obtained from the real and imaginary parts, 
respectively. 

While the general solution of the RG equation does not determine the scaling function <i> in eqn. ( V.20 ) the fact 
that it is only a function of the scaling ratio iS, allows to extract the quasiparticle structure. The analytic continuation 
s —iuj + 0+ leads to the analytic continuation of the scaling variable d —vo^ — i sign(ii7) 0+ with 



w 



(V.24) 



Writing the scaling function $ analytically continued in terms of the real and imaginary parts 4>(i? = —zu'^ — 
i sign(zi7) 0+) = ^r{w) + i ^i{w), the position of the quasiparticle pole corresponds to the value of vj for which 
the real part vanishes. Call this dimensionless real number va* , hence it is clear that the dispersion relation for the 
quasiparticles obeys 



LUr, 



VJ* v{fl) /i^ p"" 



(V.25) 



Furthermore, assuming that vanishes linearly at zu* we can write the Green's function near the position of the 
pole in the form 



Giu,p) 



1 



1 



2-r, 



(to - n7*)$^(n7*) + i^i{w*) 



(V.26) 



Alternatively, we can write the RG improved propagator, near the quasiparticle pole in the Breit-Wigner form 



GBw{i^,p) 



oj — ojp + i Tp 

with the dispersion relation, residue at the quasiparticle pole and quasiparticle width given by 

UJp — vo* v{[i) [1^^^ p^ 

Vg = {Z)W* V{fl) ^l^-^ 



Zp — 



2-)) 



(to*) 



$/(to*) w(/i) /i^"^ p'- _ Wp$/(TO*) 



to*$^(to*) 



(V.27) 



(V.28) 
(V.29) 

(V.30) 



(V.31) 



The imaginary part ^/(to) must be proportional to the anomalous dimensions, hence perturbatively small in the e 
expansion (this will be seen explicitly below to lowest order). 
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The definite values for w* , <I>i^(tn*); ^jivj*) must be found by an explicit calculation. However, the above quasipar- 
ticle properties, such as the position of the pole, group velocity, residue and width are universal in the sense that they 
only depend on the fixed point theory. For a positive dynamical exponent z the above analysis reveals a vanishing 
group velocity and width for long- wavelength quasipart icles a t the critical point. 

Furthermore the expression for the width given by (V.31) not only displays the phenomenon of critical slowing 
down, i.e, the width of the quasiparticle vanishes in the long-wavelength limit, but also the validity of the quasiparticle 
picture, since Tp/cUp <C 1 in the e expansion. 



Threshold singularities: While we have assumed above that the real part of the scaling function vanishes linearly 
at the quasiparticle pole, this need not be the general situation. It is possible that the real part vanishes with an 
anomalous power law, i.e, 



X 



X 



(2) 



(V.32) 



In this case a quasiparticle width cannot be defined as the residue will either vanish or diverge depending on the sign 
of X- It is also possible that $/(ri) also vanishes with an anomalous power at H.* . We refer to these cases as threshold 
singularities and we will find below an example of this case. Another example of this situation has been found in 
dense QCD as a result of the breakdown of the Fermi liquid theory in the normal phase|3^. Clearly only a detailed 
calculation of the scaling functions can reveal whether it vanishes linearly or with an anomalous power law at Q*. 
The set of quasiparticle properties given above ( V.28| - V.31) are only valid provided the real part vanishes linearly. 



We can go further and find the explicit form of the scaling function by foc using on the rcnormalization group 
improved propagator obtained in lowest order in the e expansion given by eqn. ( V.22 ). 

The analytic continuation to real frequencies of the RG improved two point function to lowest order given by ( V.22| ) 
leads to 



1 - 



1 + 1— sign(cj) 



(V.33) 



where we have approximated cos(^) 



sin(^) « ^ to lowest order. 



From this expression we see that th e dispersion relation ujp, group velocity Vg{p) width T{p) quasiparticles residue 
Zp are of the form given by equations ( V.2^V\3l| ) with w* — 1, and with the following explicit expressions to leading 
order in the e expansion 



z /i"'"^^ p^ 



— V = — 



(V.34) 
(V.35) 
(V.36) 

(V.37) 



with z given by ( V.23 ) above. 

Several important features of these expressions must be highlighted 



• The dispersion relation (V.25) features an anomalous dimension given by the dynamical exponent z w l+rit /2 — 
1 + e /27. The product mjT)Ji^~^ is a renormalization group invariant as can be seen from equations (V.3) and 
( |V.6| ) evaluated at the fixed point. Thus, all of the above quantities that describe the physical quasiparticle 
properties are manifestly renormalization group invariant. 



• The group velocity ( V.35| ) vanishes in the long wavelength limit as a power law completely determined by 
the dynamical anomalous dimension z. This feature highlights the collective aspects of the long- wavelength 
excitations. 

• Critical slowing down: is explicitly manifest in the width r{p) since T{p) — * as p — + 0. Furthermore we 
also emphasize the validity of the quasiparticle picture, the ratio T{p)/ujp « T^rit/2 ^ 0(e) <C 1. Thus, the 
quasiparticles are narrow in the sense that their width is much smaller than the position of the pole. Even 
considering e — 1 corresponding to dynamical critical phenomena in three spatial dimensions Tp/ujp ^ 0.1 

Thus, we see that the renormalization group resummation has led to a consistent quasiparticle picture, but in terms 
of a dispersion relation that features an anomalous dimension and a group velocity that vanishes in the long-wavelength 
limit. Obviously these features of the quasiparticles cannot be extracted from a naive perturbative expansion. 
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B. Away from the critical point: T >Tc 



Having studied the quasiparticle aspects at the critical point, we now turn our attention to their study shghtly 
away from the critical point. The critical region of interest is |T — Tel ^ ^c- Critical behavior in the broken symmetry 
phase near the critical point with T < T^, will be studied elsewhere with particular attention on the critical dynamics 
of Goldstone bosons. In this article we restrict our attention to the normal phase near the critical point. 

In the Lagrangian density (11.4) the term M^(T) is the exact mass (rather the exact inverse susceptibility) deter- 
mined by the condition (11.5) and the counterterm 5m?{T) is adjusted consistently in perturbation theory to fulfil 
this condition. However, to relate the mass to the departure away from the critical temperature it is more convenient 
to re- arrange the perturbative expansion in a manner that displays explicitly the departure from Tc- This is achieved 
as follows. Consider the one loop contribution to the self-energy in the massless theory in d -spatial dimensions: 



dfiq 



1 



Ar''^M(d) + z.T.t 



(V.38) 



with A{d) only depends on the spatial dimensionality and z.T.t stands for zero temperature terms. It is convenient 
to group this contribution with the bare mass term in the Lagrangian in the form 



= XT'^-^A{d) - TO^(O) « a{T - T^) for T^T^^ 



XA{d) 



(V.39) 



where the zero temperature contributions (denoted by z.T.t in eqn. (V.38)) had been absorbed in TO|j,(0). We 
re-organize the perturbative expansion by re-writing the Euclidean Lagrangian in the form 



1 (a.$) 



i(v$)2 



4! 2 4! 



(V.40) 



where now the counterterm 5m?' {T) is simply the one loop tadpole diagram evaluated for zero mass given by ( V.3^ ), 
it is of 0{\) and is included in the perturbative expansion consistently. 
To one-loop order the two point function is now given by 



^ +m2(T) 



I]s = (S-W(T)) = ^5: j 



d'^q 



(V.41) 



The integral above is of the typical form as those studied in the previous sections (see ( III.3| )). Separating the m = 
Matsubara contribution from the m 7^ for which we can set w?{T) cx (T — T^) « for T — Tc <C we obtain 



d'^q 



{2nY{q^+u:l^)[q^+ul + m^{T)\ 



r(^)r(i 2) ^~e(y) ^(^) rj.- 
(4^)ir(|) ^ ' ^ ' 



(V.42) 



Again, for e > we can neglect the second term in the brackets in the limit T m{T). Expanding in e to obte 
the lowest order contribution consistently in the e expansion, we obtain the two point function to one loop order 



r(2)(_ 



m^{T) 



The renormalized mass parameter mii{T) is defined by 

m^{T) = Z„i m\{T) 

and Zm is fixed by the renormalization condition 

Since Z^ receives corrections at 0{g^) we choose Zm to lowest order in the e expansion to be 



0{g\ge) 



(V.43) 

(V.44) 
(V.45) 

(V.46) 
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1. Static aspects: 

Before we embark on a full discussion of the dynamics away from the critical point, it proves convenient and 
illuminating to discuss the static aspects first. In particular since we will study the p = case but T ^ a. relevant 
quantity is the inverse susceptibility x~^{T), which is defined as 



X-\T)^M\T)^Tl{p = Q-s^Q) 



which near the critical point and the non-trivial fixed point g* given by eqn. (V.14) is given by 



K(r)] 



(V.47) 



(V.48) 



where we anticipated an exponentiation of the leading logarithms via the renormalization group, which will be borne 
out by the renormalization group analysis below. Recalhng that m?{T) cx \T ~Tc\ by eqn. ( V.39D , we find 



(V.49) 



The critical exponent 7 is seen to be the correct one Q ^ ^ p3[ . 

Just as in the case of the theory at the critical point studied above, we now study the dynamics of the theory in an 
e expansion and implement a resummation of the leading infrared divergences via the renormalization group. 

2. Dynamics away from the critical point 

As argued above the leading infrared behavior is obtained by setting the internal Matsubara frequencies to zero in 
the two loops self energy. In d = 4 — e spatial dimensions, the self-energy at two loops is 



1 



{2t:Y J {2ttY [q2 + m^{T)][k^ + m^{T)][{q + fc + p)2 + 4 + m^{T)] 



(V.50) 



The loop integrals are evaluated by introducing two Feynman parameters leading to 



l-e 



Jo JQ 

+ (1 - y)x(l - x)m'^(T) + yrn^iT) + yx% (V.51) 

It is convenient to separate the static contribution from the dynamical part by writing 

s- m{T)) = Y.^^\p- 0; m{T)) + t^^\p- s; m(T)) ; E^^) (p. ^. ^(y)) ^ j.(2) ^. ^^y)) _ 5.(2) q. ^(j.)) 

(V.52) 

The static contribution 2*^^^ (p; 0; to(T)) leads to wave function renormalization a renormalization of the mass and to 
0{e^) corrections to the static anomalous dimensions, which will be neglected to leading order in the e expansion. 
The second, dynamical contribution is obtained consistently in an e expansion: the regions of the integrals in the 
Feynman parameters that lead to inverse powers of e in an e expansion are a; ~ 0, 1 and y ~ 0. The contributions of 
these regions can be isolated by partial integration, and after some straightforward algebra we find 



S(2)(p;s;m(T)) = 



5!M 
6e 



^ -f 2^ In 



2 m^(r) -f ^ In 1 



VQm?{T) 



+ 0{e°,e) (V.53) 



Thus, putting together the one loop contribution found previously and the two loop contribution found above, the 
two point function at zero spatial momentum but away from the critical point is found to be 



1 



12e 



2 0\ 



+Oig',gh) 



2 m'{T) + 



In 1 



(V.54) 
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where we have neglected logarithmic corrections that will exponentiate to anomalous dimensions of ©(e^) for momen- 
tum dependent terms. We have only displayed the contribution that will be cancelled by wave function rcnormalization 
just as in the critical case. 



Obviously the m (T) = limit coincides with the two point function at the critical point ( [V.8 ) to leading order 



0{e). In the above expression we have not included the two loop contribution to the static s = part, since it will 
lead to an O(e^) correction to the critical exponent for the correlation length (inverse susceptibility). 

The rcnormalization conditions for the two point function away from the critical point are now summarized as 
follows: 



dT 



(2) 



= 1 



dr 



(2) 



9s2 



(V.55) 



(V.56) 



along with the rcnormalization conditions on the fou r poin t f unctio n ( [V.ll ). T o leading order in the e expansion the 
rcnormalization constants Ztj, , Zy , Zm are given by ( IV.14 ),( [V.15 ) and ( lv.46 ) respectively. 

Thus, we find the renormalized two point function at two loop order and to leading order in the e expansion (since 



In 



3e 



^UT) + -2 In 1 



vImliT) 



_9M 

o(.9^<?^e) 



mliT) 



Since m^(T) has dimension two it is convenient to introduce the dimensionless quantity 



mliT) 



and the corresponding rcnormalization group beta function 



In 



ma,T,Xo 



g + 0{g\ge) 



where we have used ( V.44 ) and ( V.46 ). 

The renormalization group equation for the A^-point function away from the critical point is now given by 



d 



d 



M-^- + Pg-^ + Pv— + (in 

an og 



d_ 

' dv 



d 
dfh? 



N 



( Si S2 SAT _ \ 

\ V V V / 



(V.57) 

(V.58) 

(V.59) 
(V.60) 

(V.61) 



The new ingredient as compared to the critical case (V.l) is the dependence on m. Following the same steps as for 
the critical case, we now find that the solution of the renormalization group equation for the two point function obeys 



r(2) 



with m(t) the solution of the differential equation 



dt 



f3Mt),v{t),fh{t)) 



with the initial condition 



-2(0) ^ rnUT (.) 



cx \T-T,{fi)\ 



(V.62) 



(V.63) 



(V.64) 
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In the infrared the couphng is driven to the non-trivial fixed point g* = e/3 and 



(V.65) 



Just as in the solution of the renormalization group equation at criticality near the fixed point ( V.18, V.19), introducing 
pe* = P ; se* = 5* we now find 



(V.66) 



Following the analysis of the critical case, and the scaling property ( |V.8| ) and writing P = /xe* we find the following 
scaling form 



S 





p' 


2-ri , 















^ is therefore identified with the correlation length^, |T[ ^ p3| 

I 



2 - 7 *-. 2 12 



(V.67) 



(V.68) 



It is important to note at this stage that the correlation length ^ is a renormalization group invariant, as can be easily 
checked by using (V.58) with the renormalization group beta function ( V.59D . 

To study the limit of zero spatial momentum it is more convenient to rewrite the above scaling solution in the 
following form 



r(2)(p,|,.g,M]=M^ 



i 



(V.69) 



From the definition of the inverse susceptibility AP{T) = x ^ (T) = F^^^ (p — 0, s = 0) we find the known result |2^, 



X-\T)^\T-T, 



7 = j—r = ^(2 - v) 



(V.70) 



Furthermore the two-point function is a function of two renormalization group invariant, dimcnsionless scaling vari- 
ables 



F(2)(p,s,m|(r,/i)) = 



to|(T, fi) 



*(^,<5) 



(V.71) 



with 



™|(T) 



2 r^2 



(V.72) 
(V.73) 



The renormalization condition (V.56) determines that \E'(0, 0) = 1. 

We can now follow the arguments provided in the previous subsection for the critical case. Under the analytic 
continuation — > —lu'^ — i sign(cL;) 0^ 



(fi -n'^ - isign(f]) 0^ 



(V.74) 
(V.75) 
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The position of the quasiparticle pole in the two point Green's function corresponds to the value of = for 
which ji{Q*{6), 6) — 0. This condition determines the dispersion relation of the quasiparticle and is given by 



(V.76) 



this expression emphasizes that the dispersion relation depends on p through the scaling variable S — (p^)^. 
Assuming that near the quasiparticle pole '^r vanishes linearly the Green's function can be approximated by 



G(p,w,mfl(T)) 



(V.77) 



where ^'ji{n*,6) = d^R{n,S)/dn\n=n' 
Breit-Wigner form 



{n~n*)^'jj{n*,5) + i'i'i{n*,5) 

Near the quasiparticle pole we can further write the above expression in the 



GBwiP,^,mii{T)) 



U! - 



i r„ 



with 



z.„ = 



cc \T-TJ 



(V.78) 

(V.79) 
(V.80) 
(V.81) 



where we have suppressed the dependence on the scaling variable S in the arguments of the real and imaginary parts 
to avoid cluttering of notation. Furthermore, we have made explicit the combination of static and dynamic critical 
exponents using the expression given in eqn. (V.68) for the static critical exponent ly and the dependence on the 
momentum is implicit through the dependence on the scaling variable S of fl* as well as the explicit dependence of 
the real and imaginary parts. 

Again, the imaginary part must be proportional t o the anomalous dimensions, hence perturbatively small in the e 
expansion. Therefore the expression for the width ( V.81 ) reveals both critical slowing down, since Tp ^ \T ~ TcY'^ 
vanishing at T = T^, and the validity of the quasiparticle picture since Tp/ujp <C 1 in the e expansion. 

At this point we recognize a fundamental difference with the Wilsonian results of ref erenc e [p^. While in ref. jl^ 
the width was found to be proportional to \T — TcY" up to logarithms, we see from (V.81) that the quasiparticle 
width actually involves the new dynamical anomalous exponent z. The difference can be traced to the fact that 
the Wilsonian approach advocated in ref.|jl^ does not include two loop diagrams which are necessary to reveal the 
anisotropic renormalization through the renormalization of the speed of light and are directly responsible for the new 
dynamical anomalous exponent z. 

We emphasize that the above Breit-Wigner form as well as the quasiparticle properties rely on the assumption that 
the real part of the scaling function vanishes linearly near the quasiparticle pole. As emphasized before in the critical 
case this need not be the general situation, and anomalous power laws can lead to threshold singularities as discussed 
above. 

While the solution of the renormalization group leads to a scaling form of the two point correlation function, it 
does not explicitly specify the scaling function ^E*. However, we can obtain the function 5* by matching the leading 
logarithms to those of the perturbative expression ( V.57) evaluated at the fixed point g* = e/3 to lowest order in the 
e expansion. Matching the leading logarithms and assuming their exponentiation via the renormalization group it is 
straightforward to see that the two point function is given by 



{s + [i + ^f-'} 



where we have used the lowest order results in the e expansion 



1 + - ; v^Oie') 



(V.82) 



(V.83) 



and kept consistently the lowest 0(e) in the exponentiation of the leading logarithms leading to ( |V.82 ). Thus, we 
obtain the lowest order result for the scaling function 



*((/3,(5) ^S + [l + (p] 



2-z 



(V.84) 
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We can now obtain an explicit form of the real and imaginary parts of the scali ng function that enter in the quasiparticle 
parameters. This is achieved by performing the analytic continuation (V.74) which leads to 



-^nin) = 6 + [l - - i sign{n) 0+] ^ 

= s-\n^ - [1 + i TT (z - 1) sign(f7)e(r22 _ i)] 



(V.85) 



For p = i.e, 5 = we see that both the real and imaginary part of the scaling function vanish at fi* = 1 with an 
anomalous power law providing an explicit example of the case of threshold singularities mentioned above. 
For p ^ and T ^ we find a quasiparticle pole at 



n* 



1 + 6— ^ l + 



2z 



(V.86) 



where we have approximated the anomalous dimension by its leading order in e using z = 1 + e/27. From this 
expression for O* we obtain the dispersion relation for quasiparticles 







2zv 






1 


[ J 


+ 




1 



(V.87) 



with V given by eqn. ( V.68), in particular we find that the frequency of zero momentum quasiparticles Wp^p |T — 
Tel^'^. Obviously at T — Tc (mij(T) = 0) the dispersion relation coincides with that of the critical case given by eqn. 
( V.25). For p 7^ 0, i.e, (5^0 the real part of the scaling function vanishes linearly and the Breit-Wigner approximation 
(V.78) near the quasiparticle pole is valid and the relati ons (V.79-V.81) describe the properties of the quasiparticles. 
To lowest order in the e expansion we find, using ( V.23) that 



2z 



(V.88) 



for p = 0, i.e, (5 = this ratio vanishes and ^'^(ri*) cx jfi* — llfj^li diverges displaying the phenomenon of threshold 
singularity with a divergent residue -Ep. 

For p ^ but T —> Tc (5 ^ oo ) this ratio equals that of the critical case (see eqn. (V.36)). 

For T ^ Tc;p ^ we finally find the width of the long wavelength quasiparticles to be given to lowest order in the 
e expansion by 



with the following behavior to lowest order in the e expansion 



p^ for p fixed, T —t 
p2z^2 ^ fixed,p^O 



(V.89) 



(V.90) 



Thus, critical slowing down emerges in both limits, furthermore the validity of the quasiparticle picture is warranted 
in the e expansion, since rjt ~ 2(z — 1) = 2e/27 + C(e^) ^ 1. 



VI. SUMMARY OF RESULTS: 



Critical phenomena, both static and dynamic in quantum field theory at finite temperature results in dimensional 
reduction since momenta and frequencies are p,uj <^ T and the correlation length is ^ 3> 1/T. The infrared physics is 
dominated by the contribution of the zero Matsubara frequency in internal loops, which in turn results in an effective 
coupling AT in the perturbative expansion. Naive perturbation theory at high temperature breaks down in four space- 
time dimensions because of the strong infrared behavior of loop diagrams near the critical point for long-wavelength 
phenomena. 

We propose an implementation of the renormalization group to study dynamical critical phenomena which hinges 
upon two main ingredients: 
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The leading infrared behavior near the critical point is determined by keeping only the zero Matsubara internal 
frequency in the loops. To control the infrared consistently we implement an expansion in e in 5 — e space-time 
dimensions. Dimensional reduction for long-wavelength phenomena near the critical point results in that the 
perturbative expansion is in terms of oc XTfi~'^ where fj. is the scale of external momenta and frequencies 
in the diagram. The renormalizcd effective coupling is driven to a fixed point in the infrared which is of 0{e). 
Therefore long- wavelength phenomena can be studied in perturbation theory around this fixed point for e <C 1. 
The perturbative expansion is improved by implementing a renormalization group resummation which reveals 
dynamical scaling phenomena with anomalous dimensions. Eventually the limit of physical interest e ^ 1 must 
be studied by further Borel and/or Pade resummations. 

The second important ingredient is the anisotropic scaling between space and time. While space is infinite, 
at finite temperature in the Euclidean formulation the time direction is compactified to the interval [0, 1/T]. 
We introduce a new parameter, the effective speed of light in the medium, which is renormalized and runs 
with the renormalization transformations. The infrared renormalization of the speed of light results in a new 
dynamical anomalous exponent which determines the dispersion relation and all the quasiparticle properties. 
The e expansion combined with the renormalization group leads to a consistent quasiparticle description of 
long- wavelength excitations near the critical point. 



The critical exponents, both static and dynamic are summarized below for the critical case T 
T ^ Tc but in the symmetric phase with T — > . 



T.. as well as for 



Table 1 : Quasiparticles at T = Tc 


ujp oc 




Vg oc 


pZ-1 


Zp oc 


pZ+ri-2 


Tp a 


Vt 


ry = 


^ + 0(e3) (static) 


m = 


^ + 0{e^) (dynamic) 


z = 


^ + H^t - v) - ^ + 4f + C(e^) (dynamic) 



(VI.l) 



Table 2 : Quasiparticles at T > Tc 


oc 




Tp oc 


'/* 1 + (p5)2^ 


Ccx 




mliT) cx 


\T~T,\ 


ly = 


^ + -^ + 0{e'^) (static) 


Vt = 


If + 0(e2) (dynamic) 


z — 


=i 1 + ^(»7t - f?) = 1 + ^ + C(e^) (dynamic) 



(VI.2) 



The new dynamical exponent z is missed by the Wilsonian approach advocated in reference ||T^ since two loops 
diagrams are completely neglected in that approach and anisotropic rescaling of frequency and momenta becomes 
manifest at two loop order and beyond. 

Critical exponents for 0{N) symmetry: 

At this stage we can generalize our results to the case of a scalar theory with 0{N) symmetry at or sli ghtl y above 
the critical point. While the static critical exponents for the 0{N) case are available in the literature|4[ pT|, p2| p3|, p5| , 
the dynamical critical exponent to lowest order in e can be obtained simply by recognizing that the symmetry factors 
correspondi ng to the 0{N) theory multiply the two loop expression for the self-energy by an overall factor. From the 
expression ( [V.g ) we see that the coefficient of s'^/v'^ is a factor 4/e times the coefficient of p^, which immediately 
leads to the result 



Vt 



4 

-V 



(VI.3) 



Since for the 0{N) theory rj = e'^ {N + 2)/[2{N + 8)^] + 0{e^) we find to lowest order in e 



Vt 



2{N + 2 

{N + sy 
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In summary the static and dynamic critical exponents to lowest order in the e expansion for the 0{N) theory are 
given by 



Table 3 : Critical exponents for 0(N). 


V — 


5 + + (static) 


V = 


2Tfw + ^(^') (static) 


Vt = 


(dynamic) 


z = 


1 + ^ WTW + ^(^') (dynamic) 



VII. CONCLUSIONS, DISCUSSION AND IMPLICATIONS 



We have studied the dynamical aspects of long- wavelength (collective) excitations at and near the critical point in 
scalar quantum field theories at high temperatures. After recognizing that naive perturbation theory breaks down 
at high temperature in the long- wavelength limit, we introduced an e expansion around 5 — e space time dimensions 
combined with the renormalization group at high temperature to resum the perturbative series. 

The effective long-wavelength theory at high temperature is described by a non-trivial fixed point at which the 
correlation functions feature scaling behavior. The anisotropy between spatial and time coordinates in Euclidean 
space-time at finite temperature leads to consider the renormalization of the speed of light, which, in turn leads to a 
new dynamical exponent z. All dynamical quantities, such as the dispersion relation and widths of long-wavelength 
quasiparticle (collective) excitations depend on this new dynamical exponent, as well as the static exponents. 

Our results are summarized in the tables in the previous section. 

Two very important aspects emerge from this treatment: i) critical slowing down, i.e, the relaxation rate of the 
quasiparticle vanishes in the long-wavelength limit or at the critical point with definite anomalous dimensions de- 
termined by the new dynamical exponent z and ii) the quasiparticle picture, i.e, narrow widths Tp <C Wp is valid. 
The group velocity of quasiparticles vanishes at the critical point in the long-wavelength limit revealing the collective 
aspects of these excitations. The dynamical exponent z = 1 + e (^^^g^i + C(e^) describes a new universality class for 
dynamical critical phenomena in quantum field theory. 

As mentioned in the introduction these phenomena have phenomenological implications for the chiral phase tran- 
sition in the Quark Gluon plasma with potential observational consequences if long-wavelength pion fluctuations 
freeze out at the chiral phase transition. An important aspect revealed by this program is that the effective coupling 
\effTn~'^ is driven to the Wilson-Fischer fixed point in the infrared, this in turn means that in this limit Ae// 0. 
This may be important in the linear Sigma model description of low energy QCD near the critical point and may give 
rise to interesting phenomenological consequences. 

In this article we focused our attention on the approach to the critical temperature from above, therefore our 
results regarding the dynamical exponent z are valid in the symmetric phase. An important question that we are 
currently addressingjS^ is the relaxation of pions slightly below Tc- Since the scattering amplitude of pions (at 
zero temperature) vanishes in the long-wavelength limit we expect novel behavior of critical slowing down for pion 
fluctuations below the critical temperature. 

We expect to report on our findings on these and other related issues soonp7|. 

While we have provided a quantitative implementation of the program of the e expansion with the resummation via 
the renormalization group, the physical limit e — > 1 requires higher order calculations with Borel or Fade resummations 
much in the same way as in static critical phenomena. We have studied the dynamical aspects to lowest order in the 
e expansion but clearly a formal proof of the consistency of the e expansion to higher orders, just as in usual critical 
phenomena, must be explored. 

While clearly such programs are beyond the scope and goals of this article, we here provided the first steps of the 
program whose potential phenomenological implications as well as intrinsic interest in finite temperature quantum 
field theory warrant further study. 
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APPENDIX A: ONE-LOOP DIAGRAM AT HIGH TEMPERATURE 



We derive in this appendix the behaviour of the one-loop diagram H{p, s) contributing to the four points function 
for high temperatures T ^ p, s. Setting w = 1 to avoid cluttering of notation (the velocity of light is not relevant for 
the discussion in this section), we have from eq.(II.7) 



Hip, s) 



T 



E 



{2T:)d- [g2 + (27r T If] [(<7 + (27r Tf {n + 0^ 



(A.Al) 



where d = 4 — e is the number of spatial dimensions. The denominators in eq. JA.AI ) can be combined using Feynman 
parameters with the result 



dx 



[q^ + Ai{x,p, s)Y 



2(4 



TT^rf2-^) f dx [Mx,p,s)]i- 



(A.A2) 



where we integrated over the spatial momenta and 

Ai{x,p,s) = x{l - x){p^ + s^) + (uji +xsf ; uji = 2ttTI] s ^ 2TrTn. (A.A3) 
We single out now the contribution from the / = mode and study the behaviour of the sum over I =/= for large 

T:^p,s. 

Let us first evaluate the I ~ term in the sum ( A.A2| ). 



T 



2(47r) 



dx [Ao{x,p,s)Y 



2(4 



^r(l) [x^Ux [(i-.y + .T^ 



(An) 



F 



2 p2 



(A.A4) 



where F(a, 6; c; z) stands for the hypergeometric function |p3|. 
We have for the I ^ terms in the high temperature limit, 

[Ai{x,p,s)]-^- ^''^■^ (2^T|/|)-^ -sign(0 



2|27rr/|2+<^ 

The sum over / ^ then yields in the high temperature limit 



{ (1 - x)p^ +s^[l-{2 + e)x]]+o{\Tl\ -3-') 



(A.A5) 



2(4 



■^r(^) f dx Y,[Ai{x,p,s)] 



3-e\ 



2 (47r)2 



12 7r(27rr)^ 



(A.A6) 



We see that only the first term in the r.h.s. is important for < e < 1 and high temperature. This term is the 
dominant high temperature limit of the sum of nonzero Matsubara modes. We then have, 



H{p,s) 



T:s>p,s 



Hasi{P,s) - 



r(l + f) C(2 + e) 



where 



Hasi{p,s) 



384 TT^'+i ri+'^ 

TAi-^r(f-i) fs^+p^ 



[p'-s\2-e)] 



l + O 



9 9 

p^, 

2^2 



A^^--^ r (1 + §) C(e) 
8 7r2+f e 



{Any 



F 



^,l-^;2. ' 



2'p2 



T 



l-e 



(A.A7) 
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APPENDIX B: TWO LOOP DIAGRAM AT HIGH TEMPERATURE 

The renormalized two points function is given by (all quantities are renormalized below) 



(B.Bl) 



where s) is the two-loops self-energy. Using the renormalization conditions ( IV.IS ) wc find for the wave 

function and velocity of light renormalizations, 



1+ V' 



0(A3) 



-fl,S—fl V 

The two loops contribution to the self energy s) is given by 



6 



(B.B2) 



(B.B3) 



where ujj = 2'k j T; s = 2t:T n . 

We combine the propagators in cq.(B.B3) using Feynman parameters and integrate over the momenta with the 
result 



\^T^ r(£- 1) 

6(47r)4-^ 



J2 dx d^ 



l-x + xai-0^ V V 



[1-x+xai-oY 



[f x^ + l' x{l-0 + U + l + nf{l-x)] 



(B.B4) 



l-e 



Using the definition ( IV. 6 ) for the dimensionless coupling, to order we find from eqs. ( B.B2 ) and ( B.B4 ) 

xHi^x)ai-0] 



= 1- 



6 \2ttT 
x{l - x)C{l - / vfJ- V 



Tie) dx di 



1-x + xCil-O 

l-2e 



Z,, — 



1 - — 

3 



5' /27rT\ 



Vfl 



[i-x+xai-or-' 
[i~x+xai-of~"' 



1 - a; -f xe(l - 



(B.B5) 



We find from the definition of the anomalous dimension (V^) and ( B.B?: ) for 7(5, — , u), 



, T ^ 2^(6+1) 



J2 dx d^ 



l-x + 2x^{l-^) 



1,3&Z 
Vfl 



xi(l-x)H{l-0 



[l-x + x^l-S)]'' 



(B.B6) 



i-x + xi{i-i) \ Vfl J 



f xC + f x(l + + l + ^) ' (1 - a^) 



-l-e 



+ Oig' 



We split the expression for 7(17, ^, v) as follows, 



T T 

7(5, = lo{9,v)+-fnz{g,—,v) 
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where 70(5,1') is the contribution from the zero Matsubara mode in eq.(B.B6) 



7o(ff,w) = g 



dx / 



[1-X+ xai - ov'^ [i~x+2 x^i - or 



and ^nz{g, ^, v) stands for the contribution of the non-zero Matsubara modes. 

For T ^ /i, we see from eq.( B.B(^ ) that jmig, ^, decreases as . [Notice that the coefficient of (j^ 

vanishes by symmetry when summing over j + I]. 

Therefore, 70(5, v) dominates for T ^ ^. 70(5, v) can be easily computed for small e > with the result. 



^2 /•! 



3+0{eg ,g ) 



Therefore, 



7o(g,w) = — / dx I d£ 

3 Jo Jo [l-x + x^{l-0] 

= ^l + 0{eg',g') 



fi i ^ 



(B.B7) 
(B.B8) 



To the lowest non-trivial order in g, that is g^ (two-loops), we find for the function /?u((7, j;;,^) 

Mg, -,v) = -| 7(3, -,v)~l - 2ej logZ. . 



(B.B9) 



where the derivatives are now at constant (bare) g. Using eq.(B.B5) for log Zy yields 



W^[,^-2e]losZy^g^^-^ 



^^(1 - 



J2 dx d^ 



X2 1(1 — x) 



[l-x + x^l-OY 



l-x + x^{1^0 



(l-^)g 



2e{l-x) 



U + i) 



T 



(B.BIO) 



where, 



Q,,{x. ^ f + (^) \fx^ + l'x{l-0]+{l- X) 

l-x + x^{l-^) \ V^J. J 



2ttT 

Vfl 



(j + + 1 



We find in the high temperature limit T ^ ^ that this expression is dominated by its zero mode contribution Wo 
(corresponding to j = Z = 0) , 



Wo = g' 



r(i + e) 



dx I d^ 



a;t-i(l - xY 



l-x + x^{l-0 



1 



lo [l^x + x^l-O] 

which turns out to be T-independent. 

The sum of non-zero terms gives a subdominant contribution for T ^ ^ and e strictly positive. We find from 
eq.(B.BlC) after calculation, 



Wn 



T»M 2 (vf^Y' r(i-i-e) 



\2TrTJ 







dx / 



X2 ^{1 — xY 



V 

,4:^. [7-2(1 



J 



ll+e 



[1-x+xai-oY 

1 + 0{^,g) . 



^^^[f(^l-^ + x^) + l2^_2jlxCY 
For < e ^ 1 and for T ^ fi, Wo and therefore W are dominated by the pole of Wo at e = 0. That is 

^T»M,0<e«l 2^2 

3e 



(B.Bll) 



(B.B12) 



where we used that x^ 



-1 e^tO 2 



Six). 



Therefore, we find for Py{g, ^,v) from eqs.( |B.B9| ), ( |B.B10D and ( |B.B12| ), 



„ , r T>^,o<£«;i V g^ 

Pv{g,—,v) = — - 



O(e0) 
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2e 
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APPENDIX C: FORMAL PROOF. 



The formal proof to one loop order begins with the expression (A.Al) from appendix |A| above. 



We now use the identity! 18, |19| 



I = T 



l — OQ 

E 

1— — 00 



[Ai{x,p,s)] 5 = 



47r i 



[A{kQ,x,p, s)] 2 coth 



fcp 

2r 



(C.Cl) 



with A{k[),x,p, s) = Ai(x,p, s;wi = —iko) and the contour C is displayed in figure below. The function 
[A{ko,x,p, s)] ^ has a cut running parallel to the real axis which for p 7^ or in the massive case for arbitrary 
p begins away from the imaginary axis and the contour C. The contour can now be deformed to wrap around the cut 
and the analytic continuation s — > —iuj + 0"*" can be performed. For p,uj -^T the infrared behavior is dominated by 
ko ^ T for which coth[A;o/2T] ~ 2T/ko and the resulting expression features a pole at fco = while the cut begins 
away from the origin. The cut can be deformed again to circle the origin and the integral is simply the residue at 
the pole fco = 0. Therefore the infrared dominant term is given by lir = [Ao{x,v,s = —ioj + O^)] ^, result which 
coincides with the analysis in terms of the Matsubara sums provided in appendix [AJ. 



[ko 



FIG. 1: Contour in the complex fco plane 
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